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TECHNICAL PUBLICATION 


OPTIMIZATION OF CREW SHIELDING REQUIREMENT 
IN REACTOR-POWERED LUNAR SURFACE MISSIONS 

1. INTRODUCTION 


In addition to other flight risks and hazards, space flight beyond the confines of the Earth’s mag- 
netic field must face the challenges of space radiation exposure. In extended lunar surface missions, pro- 
tection of crew and systems requires shielding strategies against various sources of space radiation fields, 
both natural and man-introduced. Due to various degrees of variability, unpredictability, and, in some 
critical areas, lack of basic data, guaranteeing safe levels of exposure poses a special challenge. 

Exposure estimates for shielding solutions as well as for safety assessment must be formulated 
and optimized based on incomplete data, constrained by both technical and nontechnical factors. One 
of the more consequential constraints, albeit somewhat subjective, is that of “as low as reasonably 
achievable” (ALARA). 

A main task of mission designers is to minimize requirements on structure and function while 
ensuring maximum protection for crew and systems, consistent with ALARA. ALARA is currently 
NASA’s accepted guideline as well as being a part of the legal requirements with regard to ionizing radia- 
tion exposure and crew health and protection. 

To be consistent with ALARA, shielding solutions and dose and risk assessments must rely 
on robust and accurate exposure estimates. Objective comparisons among these solutions will clearly 
require reliable estimates as well. 

To various degrees, such estimates are hampered by inherent uncertainties in basic knowledge 
of the radiation environment itself, its transport and interaction in various media of complex geometry and 
composition, and most critically, in the human biological response to such exposure. 

In the absence of more empirical data, on the one hand, and the increasing complexity of the 
modality and applications by which (and for which) these estimates are determined, on the other, such 
estimates are best viewed as guidelines rather then predictions. 

Given the expected doses, this parametric study focuses on estimating the optimal crew shielding 
requirement in lunar surface missions with a nuclear option. Possible missions are assumed to take place 
during both low and high solar activity. Specificity due to the mission’s location on the lunar surface is not 
taken into account. For this study’s purposes, these missions are assumed to only include a crew habita- 
tion module that is powered by a small fission reactor placed at some distance from this module. No other 


details about the reactor or the habitation module, such as their geometric configurations and specific 
structures or subsystem, are either assumed or used. 

Independent of the exact type or chemical composition of the shielding material, any shielding 
solution will require a certain amount of areal density to reduce the crew’s expected exposure to accept- 
able levels. For this study, lunar regolith, albeit in an idealized form, is assumed to be the shielding mate- 
rial of choice . 1 

The estimates and method presented here are meant to help mission designers put in perspec- 
tive the expected cumulative exposure — due to natural and introduced sources — vis-a-vis the amount of 
regolith mass required for crew protection. For example, for logistical considerations, it may be desirable 
to minimize the separation distance between habitat and reactor while maintaining maximum protection. 
Conversely, it may be desirable to minimize the amount of regolith to be used by maximizing the distance. 
Ideally, in both extremes, as well as for all estimates in between, required regolith mass must be optimized 
for each separation distance. 

Since shielding will be required that can be used for both reactor and crew, a self-consistent 
approach would be to estimate, at a given distance, the optimal and also total amount of regolith mass 
separating the crew from the reactor. Because of the additive nature of the solution, this amount can be 
thought of as the sum of habitat and reactor shielding. This self-consistent solution should allow for more 
flexibility in allocating material resource and/or construction effort between reactor and habitat. 

A brief survey of the radiation environment and exposure doses is presented, followed by a descrip- 
tion of the dose-depth relations used and the one-dimensional optimization scheme. Sample results for 
optimized required regolith mass and reactor/crew separation for missions during solar minimum and 
maximum conditions, superimposed on a ‘typical’ large solar particle event, are provided, followed by a 
discussion and conclusions. 
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2. THE RADIATION ENVIRONMENT 


Energetic, high-charge galactic cosmic-ray (GCR) ions and solar energetic particles (SEP) consti- 
tute the main (natural) source of this intense radiation environment. The energy range of these particles 
spans more than eight orders of magnitude (from thermal to ultra-relativistic) while their atomic numbers 
populate all of the stable nuclides of the periodic table. 

Atomic charges of 1 (hydrogen) through 26 (iron), however, are considered important for crew 
radiation safety and shielding purposes. By number, hydrogen constitutes about 90%, helium, 7%, 
and all others, 3% of the GCR ions. The intensity of the ambient GCR component (=1 cm~ 2 ) peaks around 
500 MeV/nucleon and is modulated by a factor of about three over the 1 1 -year solar cycle. 2 During solar 
maximum, and due to the actions of the solar wind, access to the heliosphere by diffusing GCR ions is 
reduced. As a result, the GCR component appears depressed in the inner heliosphere. 

During heightened solar activities, solar particle events (SPE), while random in occurrence, are 
more frequent and are strong enough to transport SEPs, by a propagating shock driven by a coronal mass 
ejection (CME), to Earth’s orbit and beyond. The SEP component is mostly composed of energetic pro- 
tons, and it peaks around a few tens of MeV in energy, but can vary widely in intensity (~l O 7 cm~ 2 ) as well 
as in the shape of its energy spectra. The so-called ‘large’ events, such as the October 1989 event, can be 
an order of magnitude more intense than the ‘average’ event, and many orders of magnitude above the 
quiescent conditions, lasting for hours to 2-3 days. 2 

Relatively little is known (or can reliably be predicted) about the photospheric, coronal, and helio- 
spheric mechanisms responsible for CMEs and large SPEs. Furthering our basic understanding in these 
areas remains a key prerequisite 3 of the Exploration Vision. 

In addition to these natural sources of space energetic particles, there are likely to be man-intro- 
duced radioactive and fission sources for power and even propulsion purposes as well. A number of stud- 
ies 4 ’ 5 conducted for the power requirements during future lunar surface missions, for example, suggest 
that the need is on the order of tens of kilowatts of electric power. 

For this level of power, chemical, solar, and radioisotope sources may be insufficient or impracti- 
cal. For crew protection purposes, fission reactors are considered mainly as sources of energetic neutrons 
and gamma rays (photons). Contributions of these sources to the total expected crew dose are mostly due 
to prompt neutrons. Prompt neutrons are produced in the fission process of the fissile material, such as 
uranium-235 or plutoium-239. Most of these are energetic or ‘fast’ neutrons produced (at =10 14 cm~ 2 ) 
as direct fission products with an average energy of about 2 MeV. Photons (at =10 10 cm~ 2 ) are produced 
as both direct products of the fission reaction as well as a result of the subsequent decay of the fission 
radioactive products. For shielding purposes, however, gamma rays with energy <0.6 MeV are typically 
ignored. 6 
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3. EXPECTED EXPOSURE LEVELS 


Crew exposure levels are typically expressed in dose equivalent units. Dose equivalent, in Sievert 
(Sv) units, is calculated from the dose corrected by a dimensionless, multiplicative factor called the radia- 
tion ‘quality factor’, or Q-factor. 7 Ionizing radiation like energetic heavy ions (such as GCR ions) are 
characterized by high Q values. Uncharged neutrons are also assigned high Q values to underscore their 
more serious health hazards relative to either x rays or gamma rays at the same energy. Unlike the phys- 
ically describable and measurable dose, the Q-factor is an empirical, dimensionless variable assumed 
to ‘represent’ the majority of the biological effects associated with exposure to ionizing radiation, 
but without specifying such effects by their end points or response functions. 8 

Estimating the health risk — and thus, shielding requirement — associated with space radiation 
exposure is hampered mostly by uncertainties in the biological response. 9 Other factors associated with 
the radiation environment and its physical interactions, as well as the dose and dose-rate volatilities, 
also contribute to the risk. As will be discussed later on, large (-200%) uncertainties in the Q-factor can 
significantly affect shielding requirements, and any optimized estimates of which, as well. 

The National Commission on Radiation Protection (NCRP) publishes and regularly updates 
recommended limits appropriate for low Earth orbits (LEO) missions. Table 1 lists the 1999 recommenda- 
tions 10 for organ dose limits for all ages for 30-day, annual, and career exposures. (Note the 50-cSv limit 
for bone marrow.) 


Table 1 . 1999 NCRP recommended dose limits by organ and exposure duration. 


Limit 

(cSv) 

Bone 

Marrow 

Eye 

Skin 

30-day exposure 

25 

100 

150 

Annual 

50 

200 

300 

Career 

50-300 

400 

600 


To put this 50-cSv limit and the other NCRP limits in perspective, on the International Space 
Station (ISS), for example, during solar maximum, the average effective dose was measured to be about 
6.1 cSv, while the effective dose rate was about 0.037 cSv/day. 11 Note though that on the ISS, in addition 
to protective geomagnetic effects (which are not present outside the magnetosphere), shielding equivalent 
to about 5-10 cm of aluminum is provided by the ISS structure and system’s materials. 12 

On the lunar surface, the dose due to the (isotropic) GCR source is reduced by a half due to the 
shadow shielding effect of the Moon itself. The introduction of a small nuclear fission reactor (— 25-kWe) 
is estimated 13 to add about 5 cSv/year at a ‘safe distance’ from its shielded core. Both water and regolith 
have been considered for core shielding. 13 
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Table 2 contrasts typical expected 14 ’ 15 doses on the surface of the Moon with and without a 50-cm 
thick shield made of idealized lunar regolith, equivalent to 1 1 in of standard aluminum, assuming solar 
minimum GCR conditions and superimposed on an August 1972 class SPE. Given currently accepted limits 
for LEO missions (see table 1) these expected exposure figures clearly suggest that extended (>6 mo) 
surface missions will require shielding solutions, even without the presence of a nuclear fission source. 


Table 2. Expected doses on the lunar surface with and without 
shielding (no nuclear power source assumed). 


Duration 

GCR 

SEP 

Mission 

(days) 

(cSv) 

(cSv) 

(cSv) 

10 

0.3/0. 8 

7.5/20.5 

7.8/21.3 

30 

1/2.5 

7.5/20.5 

8.5/23 

180 

6/15 

7.5/20.5 

13.5/35.5 

360 

12/30 

7.5/20.5 

19.5/50.5 
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4. PARAMETERIZING THE DOSE-DEPTH RELATIONS 


For the purpose of this parametric study, dose as a function of depth in lunar regolith from 
all three radiation sources, i.e., GCR, SEP, and fission sources (we ignore contribution from neutron 
albedo) will be assumed to have simple closed form expressions amenable to variational analysis. 
To that end, the GCR dose-depth relation is taken to be 

D, (x) = A, exp(-A,x) + B ] , ^ 

where D | (x) is the dose rate in cSv/yr, x is total regolith separation mass between reactor and crew 
in g/cm 2 (i.e., an arbitrary combination of reactor depth and habitat shielding) and Aj is the regolith 
attenuation coefficient for GCR in (g/cm 2 ) -1 . 

The constants A 1 = 74 cSv/yr and B 1 = 28 cSv/yr, as well as A 1 = 0.06 (g/cm 2 ) -1 are estimated 
using fits to three-dimensional Monte-Carlo simulations assuming solar minimum conditions. 16 For 
solar maximum conditions, the values are: A 1 = 54, By = —24, and A 1 = 0.02. 

For this approximation as well as for the other two below, lunar regolith is idealized as being 
composed of 74% oxygen, 11% silicon, 7% aluminum, 4% calcium, and 4% magnesium by weight. 
The density of this aggregate is taken to be 1.5 g/cm 3 . 

The GCR/SEP particle flux is transported through a thick slab of this idealized regolith, suf- 
fering both energy and charge losses. The transported flux is converted into dose and dose equivalent 
quantities using the International Commission on Radiological Protection (ICRP) 1991 conversion 
convention. 17 


The SEP transported flux is similarly assumed to be of a simple (analytic) form, 


D 2 (x) = 


An 


B 2 + A 2 x 


( 2 ) 


where D 2 (x) is now the event integrated dose in cSv, A 2 = 400 cSv, B 2 = 1, and X 2 = 1.08 (g/cm 2 ) -1 . 
These numbers are based on three-dimensional transport simulations through a finite slab of lunar 
regolith as described above and for an assumed August 1972 class SPE. 

The dose-depth approximation as a function of radial distance from the reactor’s location is 
also based on three-dimensional transport simulations. 13 The conceptualized reactor in the simu- 
lation is a moderated spectrum, NaK cooled, Hastelloy™ (a Hayes International, Inc. product)/ 
uranium-zirconium hybrid (UZrH) reactor with open lattice pin geometry. 18 The reactor provides 
thermal power to a 25-kWe Stirling engine power conversion system. The cylinder- shaped system 
(reactor, water shield, and power conversion system) stands ~ 2 m high and is ~ 1 m in diameter. 
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The reactor’s transported 13 neutron and gamma rays fluxes are assumed to originate from a shielded 
core. To first order, the reactor’s dose-depth relation for a given r (surface separation distance in m) can be 
approximated as 


D 3 (x) = (A 3 exp(-X 3 x) + B 3 )/r 2 , (3) 

where D 3 (x) is the dose rate in cSv/yr, A 3 = 2 x 10 6 cSv/yr-m 2 , B 3 = 3xl0 3 cSv/yr-m 2 , and A 3 = 1 ,87x 1 0 2 
(g/cm 2 ) -1 . 
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5. OPTIMIZATION SCHEME 


To formulate a one-dimensional variational scheme, equation (3) is re-expressed as a controllable, 
‘dynamical’ system as: 


df) 3 , 5 3 

— = -fl 3 (r) + -yc(i) , (4) 

where x = A 3 x is the ‘dynamical’ variable, c(x') is the control variable, and r is a parameter. The con- 
trollability of the process is assumed 19 based on the system being autonomous, linear, and possessing 
a stable, (uncontrolled) ‘equilibrium’ state asx'— > «, 

The initial condition, D 3 (0), is taken to be the uncontrolled state atx'= x = 0 where the control vari- 
able c is identically equal to unity. Formulated this way, the objective becomes to find the optimal regolith 
mass, x'= x*, such that for a given r the functional: 


J 


[*>)]= J, 


x (r) 
0 


T + — C 2 (x / ) 
2 


dx 


(5) 


is minimal while assuring a safe dose, i.e., D 3 (x*) < D s . 

An optimal solution is assumed to exist due to the convexity property of J [x*(r)J, i.e., over its 
entire domain V, J(pc') assumes a minimum value at each and every stationary point in D. This property 
of J assures 20 that 


J(x') > J(x*) + V/(x*) ■ (x' - x*); Vx',x* e V , ^ 

where V.7 is the gradient of J. 

The safe dose, D is taken to be the dose limit (an NCRP limit) including the contributions due to 
GCR and SEP exposure as a function of depth x . The first term in this ‘cost’ functional J is taken to be 
solely determined by the total mass required, x*, while the second term by the incremental amount of mass 
needed to reduce the incurred dose to its current level at this x\r) point. 

The constant x is a measure of this distribution between the two: When x » 1 , this corresponds 
to a solution for achieving a safe dose level at a given r with as little regulation, i.e., /--manipulation, as 
possible. Conversely, when x <sc 1 , the safe dose level is achieved for maximal manipulation (regulation). 
Note that no optimal solution exists when x is identically zero. 


8 


The optimization proceeds by assigning a ‘Hamiltonian’ to the process according to the Pontryagin 
maximal principle . 19 ’ 21 The Hamiltonian remains constant along an optimal trajectory, x'= 0 — > to x = x* . 
The general form for a one-dimensional Hamiltonian is: 


H = Vq/1q + V ] fl ] . (7) 

The pi variables are called state variables while the v ones are called the costate variables (analo- 
gous to generalized coordinates and generalized momenta in analytical dynamics). Both sets are given 
by Hamilton equations of motion, 



dH 

dH 



O 

& 

i 

ii 

< 

ll 

i 

(8) 


dH 

dH 



Bo ~ + -j 

dv 0 

Bl ~ + 3 
dv. 

(9) 

At each point along the optimal trajectory the Hamiltonian remains 

minimized. For this system , 21 

i.e., equations (4) and (5), the Hamiltonian is: 



H(x') = - 

T + — C 2 (x') 

+ V\[-D 3 (x') + c(x')] 

(10) 


Solving for Vj and c, and applying initial and safety conditions onD 3 , gives the following transcen- 
dental relation for x*(r): 


£> 3 ( 0 )/^. exp x* -D s (x*)/D r - 2/3 sinh x* = 0 , (11) 

where D r = B^/r 2 and /3 is a constant of the ‘motion. ’ Constants of the motion a and /3 are determined from 
the initial and safety conditions, 


P - Z) 3 (0)-a/2 , 


( 12 ) 


and where a is the negative root of: 


a 2 -2aD 3 (0)-2r = 0 . (13) 

The next step is to estimate the value of the constant x for this particular optimization, equation (4), 
and the choice for the functional form of J, equation (5). 

The directional derivative of Jat c is defined as 20 : 


8J(c;d ) - lim 
£ — >0 


J(c + ed ) - J(c ) 
e 


( 14 ) 
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(15) 


— ~ — (C + £(^) g_Q . 


From equation (5), 


J(c + £ci) — J* 



— d 2 dx . 
2 


(16) 


Subtracting J(c), dividing by £, and taking the limit as £ — » 0, results in 


8J(c;d)= [ c(x')d(x')dx' . 

•'O 


(17) 


Now, from the convexity property of J, 



(18) 


and from the symmetry property of equation (17) with respect to c <-> d, and recalling that c = 1 corre- 
sponds to the uncontrolled, initial condition, results in 


From the general solution of equations (4-11), it is known that c(x') <=c cxp(x'). It follows then, from the 
above relation, that, to within a constant of order unity, the numerical value of T should be -By 

It should be noted that for this particular optimization scheme of equation (3), a different approach 
would have been to use the conditions on the Hamiltonian, i.e., minimal (including zero) and unchanged, 
along an optimal trajectory, rather than minimizing the cost functional J, as was done here. The alternate 
approach should, in principle, give the same results, but no attempt, for self-consistency, has been made 
here to demonstrate as much. 

The first order, linear optimization scheme presented here should also be treated as parametrization 
specific in so far as the fonn of equation (3) is concerned, i.e., its x and /'-dependence and the treatment, 
for purposes of estimating the optimal path, x — > x*, of the variable x as the ‘dynamical’ variable and r 
as being part of the control variable c{x'). No attempt has been made here to check for the applicability 
of the solution (controllability, existence, uniqueness, etc.) over wide ranges of the fit parameters, A t , B h 
and X- v However, the theory of linear, first order control problems, such as the one described by equa- 
tion (4), is well anchored, and properties of the general solutions are known for sufficiently large phase 
and parameter spaces, especially so for autonomous, one-dimensional systems. 



(19) 
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The choice of the cost functional, equation (5), also affects the solution; convexity-wise only 
the simplest form of the functional, i.e., quadratic, was used. Clearly, and as is discussed below, other 
forms must be explored as well. Finally, generalization of linear-state control problems to two and three 
dimensions is, in principle, straightforward. However, issues related to uniqueness and stability of the 
controlled solution become more critical in higher dimensions. 21 Generalization of this particular optimi- 
zation scheme to higher dimensions must be preceded by further numerical and analytical demonstrations 
of its wider applicability and utility. 
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6. SAMPLE CALCULATION AND DISCUSSION 


The optimization scheme demonstrated above is applied to two mission scenarios: one during 
GCR solar-maximum conditions superimposed on an August 1972 class SPE (figure 1), and the second 
is done similarly for GCR solar-minimum conditions (figure 2). (Note that ‘depth’ in the figures refers 
to the total regolith mass, in g/cm 2 , separating reactor from crew.) 

In both cases, the shielding material is the idealized lunar regolith as described in section 4, along 
with the parameterized forms and values of the transported radiation sources for each scenario. For each 
scenario, equations (4-11) are solved using the fit (A^ B p and A ; ) and optimization (a-, /3 ; , and T ; ) param- 
eters, self-consistently. These latter ones depend sensitively on initial conditions and hence they change 
from one scenario to the other. The dose limit, for reference, is taken to be the 50-cSv/yr level, i.e., 
the LEO 1999 NCRP annual limit for bone marrow exposure (table 1). 

For each scenario, as a function of distance from the reactor, shown in the figures is the optimized 
total (due to reactor plus natural) mass of lunar regolith required to keep the dose rate level less than 
or equal to the safe rate of 50 cSv/yr. Also shown is the required mass for the reactor-only case, i.e., no 
GCR or SEP fields assumed, and for the natural environment-only case, i.e., no reactor. After subtracting 
the mass requirement to shield against the GCR and SEP fields, the balance can, as mentioned earlier, be 
treated as an arbitrary combination of both the amount of shielding required for the reactor plus that for 
added shielding, due to the introduction of the reactor, for the habitat. 
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Figure 1 . Optimized reactor depth-distance solution for GCR solar maximum conditions 

superimposed on an August 1972 class SPE, for a 50-cSv/yr dose limit (see table 1). 



Figure 2. Optimized reactor depth-distance solution for GCR solar minimum conditions 

superimposed on an August 1972 class SPE, for a 50-cSv/yr dose limit (see table 1). 


For example, for a surface mission during solar minimum, at a distance of 100 m from the reactor, 
from figure 2, the optimized total regolith shielding requirement is about 62 g/cm 2 . Shielding against GCR 
and SEP fields requires about 16 g/cm 2 . Note that the un-optimized reactor requirement (which is also 
the total here because it is larger than the natural overburden) is about 76 g/cm 2 , which is a 23% savings 
in required mass due only to optimization. (For solar maximum conditions, figure 1, the savings are, 
of course, even larger (30-35%) because the natural environment overburden is lower.) 

In addition, the 46-g/cm 2 requirement can be divided in a number of ways depending on other 
factors, such as availability and processing of regolith and reactor site preparation, between the actual 
required depth of the reactor system beneath the lunar surface and the actual thickness of the added habitat 
protection against the reactor’s radiation fields. This added flexibility is a result of treating the reactor 
and habitat shielding requirements self-consistently in this simple optimization scheme. 

However, this self-consistent treatment is also reflected in the optimization cost. In figure 2, 
for example, and for distances larger than about 133 m from the reactor, the ‘optimized’ mass is larger 
than what is actually required. The reason being the ‘cost’ of optimizing the mass for any distance 
is always nonzero, as can be seen from equation (5). In this particular optimization scheme, the optimiza- 
tion becomes ‘cost-ineffective’ for large distances, but not large enough, i.e., for distances at which the 
reactor’s fields become negligible compared to the natural overburden (=220 m for this study). Clearly, 
a more robust fonn for the cost functional, equation (5), is required to reduce the cost over a wider range 
of separation distance. 

Also, the above assessment was based on an idealized regolith and its simulated attenuation prop- 
erties against both natural and fission radiation sources. If one allows for an error margin of the same order 
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in the attenuation properties of regolith (and not in its other physical properties 22 ), this savings all but 
disappears. Imprecision in basic regolith attenuation properties that is on the order of 50-75% will render 
any optimization scheme frivolous. 

It is important to note that variations in regolith density alone, which has a range of 1 .5-2.8 g/cm 3 , 
can easily contribute to this level of imprecision. When coupled with uncertainties in modeling the radia- 
tion quality factor, it becomes clear that this and similar optimization schemes are easily defeated by such 
large variabilities. Unfortunately, some of these variabilities are inherent to shielding and radiation protec- 
tion studies associated with crewed lunar surface missions, with or without a nuclear option. 
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7. CONCLUSIONS 


A parametric study was conducted to afford mission designers first order estimates for the amount 
of lunar regolith required to protect the crew on a lunar surface mission from exposure to GCR, SEP, 
and neutron fields associated with a small fission reactor. 

Since shielding is expected to be required for both reactor and crew, a self-consistent approach was 
taken to estimate, at a given distance, the optimal (total) amount of regolith separating crew from reactor. 
The additive nature of the solution in this treatment should allow for some flexibility in allocating material 
resource and/or construction effort between reactor and habitat. 

Simple but simulation based dose-depth relations were used for all three radiation sources in a one- 
dimensional optimization scheme. The objective was to estimate the optimal regolith mass between crew 
and reactor, as a function of their separation distance. The optimization scheme was based on Pontryagin 
maximal principle. 

The scheme was applied to both solar maximum and minimum conditions. Depending on the 
mission’s time profile, a savings of up to 30% in mass can be realized between optimized and un-optimzed 
required regolith mass estimates. However, it is argued that variation and uncertainty mainly in lunar rego- 
lith attenuation properties and in the radiation quality factor can easily defeat this and any other similar 
optimization scheme. 
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